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Abstract 

Using the Feynman-Vdovichenko combinatorial approach to the two dimensional 
Ising model, we determine the exact Curie temperature for all two dimensional 
Archimedean lattices. By means of duality, we extend our results to cover all two 
dimensional Laves lattices. For those lattices where the exact critical temperatures 
are not exactly known yet, we compare them with Monte Carlo simulations. 

1 Introduction 

As a solvable example of a statistical system exhibiting a second order phase transition, 
the Ising model has been of central importance in the development of statistical mechanics 
and still continues to play a key role as test ground for new theoretical and computational 
methods. 

The existence of a Curie point, at which a transition from the paramagnetic phase to 
the ferromagnetic phase takes place, as the temperature is lowered, was first proven to 
exist by Peierls [I] in the case of the square lattice Ising model. Successively Kramers and 
Wannier [2] calculated the value of the critical temperature using duality and the property 
of the square lattice of being self-dual. With the first calculation of the zero field free 
energy by Onsager [3J H] and the calculation of the spontaneous magnetization by Yang 
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[5], the second order character of the phase transition was established. This successful 
attempt to solve the Ising model on a square lattice was carried over using transfer matrix 
techniques. Along the same path, Wannier was first to study the triangular and hexagonal 
(also called honeycomb) lattice Ising model, in particular in relation to antiferromagnetism 
studies [6]. The first non regular lattice zero field free energy was calculated for the 
"Kagome" lattice by Kano and Naya [7J. Successively the "extended Kagome" lattice was 
studied by Soyzi [8J. The "Union Jack" lattice and the dual, the "Bathroom-tile" lattice, 
were covered by [HI HP] . Other lattices were investigated in [TT| IT2] . 

A combinatorial approach to solve the Ising model was proposed by Kac and Ward 
[14]. It is based on the possibility to map the high temperature expansion of the parti- 
tion function onto a particular random-walk problem on the lattice. Feynman [15] and 
Vodvicenko [16J showed how to recover the Onsager solution by this method. The above 
mapping is based on a topological theorem about loops on two dimensional lattices proven 
by Sherman |17| . 

Despite of many efforts, the free energy in presence of a non-zero external magnetic 
field has not yet been calculated exactly. Some results are available only at criticality 
T = T c or for weak magnetic field B — > [23]. Still nowadays the Ising model is an active 
area of research. 

In this prospect any exact analytical result can be useful. In this paper we use the 
combinatorial method to calculate the Curie temperatures for the Ising model on all 
Archimedean lattices and all their duals, the Laves lattices. This method was first used 
by Thompson and Wardrop [19J to calculate the critical temperatures for two Archimedean 
lattices, here we calculate the Curie temperature for all of them, including three values 
which were known by now only approximately from Monte Carlo simulations. 

2 Archimedean and Laves lattices 

Archimedean lattices are uniform tilings of the plane in which all the faces are regular 
polygons and the symmetry group acts transitively on the vertices. It follows that all 
vertices are equivalent and have the same coordination number z. Also, all vertices are 
shared by the same set of polygons and thus we can associate to each Archimedean lattice 
a set of integers (pi,p 2 ,...) indicating, in cyclic order, the polygons meeting at a given 
vertex [20]. When a polygon appears more than one time consecutively, for example as 
in (...,p,p, ...), we abbreviate the notation writing (...,p 2 , ...). In this way, the uniform 
tilings of the plane, made with squares, triangles and hexagons, are indicated respectively 
by (4 4 ), (3 6 ) and (6 3 ). 

It is possible to show, that apart the regular tilings just listed, there are eight other 
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(3 2 ,4,3,4) (3,4,6,4) (3,6,3,6) (3, 12 2 ) 




(4,8 2 ) (3 3 ,4 2 ) (4,6,12) (3 4 , 6) 



Figure 1: The eight uniform tilings of the plane made with more than one regular polygon. 

uniform tilings of the plane. This are all constructed using more than one regular polygon 
and are shown in Figure 1. In total we have eleven different two dimensional Archimedean 
lattices. Some of this lattices have names: "Kagome" (3,6,3,6), "extended Kagome" 
(3, 12 2 ) and "Bathroom-tile" (4, 8 2 ), while the others have no particular name. Note that 
the lattice (3 4 , 6) is the only non-specular one. 

The dual lattice A* of a given lattice A is constructed associating a vertex of A* to 
each face of A. The duals of the Archimedean lattices are called Laves lattices. They 
are tilings of the plane made with one non-regular polygon and are indicated by the 
notation [pi,p2, •••]• The square lattice is self-dual [4 4 ] = (4 4 ) (this property was used by 
Kramers and Wannier to first calculate the critical temperature for this lattice), while the 
triangular and hexagonal lattices are dual of each other [3 6 ] = (6 3 ), [6 3 ] = (3 6 ). Of the 
duals of the other uniform lattices, three are triangular tilings, two are rhomboidal tilings 
and three are pentagonal tilings. The dual of the Kagome lattice [3, 6, 3, 6] is called "Dice 
lattice" and the lattice [4,8 2 ], dual to the "Bathroom tile" lattice, is called "Union Jack" 
lattice. 
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3 Ising model 



We briefly review the definition of the model and the basic steps involved in the combi- 
natorial solution of it. At every lattice site % a spin variable oi G { — 1,1} is located and 
a microstate is defined by a spin configuration {cr}. Nearest neighbor interactions are 
assumed so that the energy of a given spin configuration is given by 

%} = -J^^-B^a, (1) 

If J > the interaction is ferromagnetic, while it is antiferromagnetic if J < 0. For the 
aim of calculating the Curie temperature, we can set to zero the external applied magnetic 
field B = 0. The partition function is the sum over all spin configurations, weighted by 
the Boltzman-Gibbs factor 

Z A = ]Texp ji^^ l . (2) 

We defined K = J/ksT where ks is the Boltzmann constant. 

Using the high temperature expansion, the partition function ([2]), can be rewritten as: 

Z A = 2 N (coshK) NE $ A (v). (3) 

In ([2]) we defined v = tannic, while iV is the number of lattice sites and N& = ~iV is 
the number of links per lattice site. The function Q\(v ) is the generating function of the 
numbers G n (we define Go = 1), which count the graphs of length n with even vertices 
that can be drawn on the lattice A. In this way, the problem of solving the Ising model 
on A is reduced to the combinatorial problem of counting even closed graphs on A. 

If we don't care about the fact that we have to count only over closed graphs with even 
vertices, we can use the standard combinatorial property that the sum over all (possibly 
disconnected) closed graphs, is given by the exponential of the sum over only connected 
closed graphs. But connected closed graphs can be seen as closed random-walk paths. 
If we find a way to consider in the sum only closed paths with even vertices, we have a 
way to calculate the generating function Q\(v) simply by counting closed random- walks 
paths. 

Feynman [15] and Vodvicenko [T6] introduced a trick to count only over paths with 
even vertices: weight the random-walk with a complex amplitude a = e td l 2 at every turn 
of angle 9 the random-walker makes along his trajectory. In this way all contributions 
coming from non-even paths delete each other in the sum, and the mapping from the 
high temperature expansion to the random-walk problem works out correctly (see [251 121] 
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A z p m 

(3,12 2 ) 3 6 18~ 

(4,6,12) 3 12 36 

(4,8 2 ) 3 4 12 

(6 3 ) 3 2 6 

(3,4,6,4) 4 6 24 

(4 4 ) 4 1 4 

(3,6,3,6) 4 3 12 

(3 4 ,6) 5 6 30 

(3 3 ,4 2 ) 5 2 10 

(3 2 ,4,3,4) 5 4 20 

(3 6 ) 6 1 6 

Table 1: Coordination number z, number of lattice sites per unit cell p and size m of the 
random-walk matrix Wa for each Archimedean lattice. 

for more details). This results have a rigorous base on a topological theorem proven by 
Sherman [T7] and were further investigated by Morita |18j . 

After introducing the appropriate amplitudes, which depend on the properties of the 
lattice A, we need to count all connected closed graphs of length n that can be draw on 
A. This number is equal to the number of random-walks of length n starting from all the 
lattice sites 

— X) iVto(xo,xo,n), (4) 

xo.io 

where iVj(x, Xq, n) is the number of walks arriving at site x from direction i at the n-th 
step starting from site xo. The sum is over all possible starting points and the factor j- 
exclude the over-counting due to the fact that there are two directions in which a graph 
can be travelled and the fact that each graph can be seen as a random-walk starting from 
all of the sites of the lattice reached by the random- walker. 

The random- walk problem can be written in terms of the master equation 

A^(x, xo, n) = W itj {x, x^A^x', x , n - 1) , (5) 

with boundary conditions A^(x, x , 0) = 5 x ,x ^i,io- Equation flS]) is solved by A^(x, x , n) = 
Wj™ (x, xo). From equation (j3J) we find, in the limit — > oo, the explicit form for the 
generating function: 

{ N f d 2 k 1 
$ A (w) = exp | — J — — log det [1 - vW A (k x , k y )} \ , 

where the integration is over the region < k x < 2tt and < k y < 2tt. The free energy 
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z 


A 


Pa(v) 


3 


(3,12 2 ) 


(1 + v) A (l - 2v + 3v 2 - 2v 3 - 2v A ) 2 




(4,6,12) 


(1 + 2v 2 + 5f 4 ) 2 (1 - 2v 2 + 2v A - 10v 6 + v 8 ) 2 




(4,8 2 ) 


(1 - d - w 2 )(l + u + 2v 2 )(l - 4v 3 - v A ) 




(6 3 ) 


(1 - 3v 2 ) 2 


4 


(3,4,6,4) 


(1 - v 2 ) 2 (l + w 2 ) 6 (l - Av 2 - 6f 4 - 4t> 6 + v 8 ) 




(3,6,3,6) 


(1 - v 2 ) 2 {\ - Av 2 - 6v A - 4v 6 + V s ) 




(4 4 ) 


(1 - 2v - V 2 ) 2 


5 


(3 4 ,6) 


(1 + v) 8 (l + 3w 2 ) 2 (l - 4v + 7w 2 - 12w 3 + 3w 4 - 3u 6 ) 2 




(3 3 ,4 2 ) 


(l + v) 2 (l -3v) 2 (l + v 2 ) 2 




(3 2 ,4,3,4) 


(1 + v) A (l -2v-v 2 - Av 3 - 9w 4 + 6v 5 - 7v 6 ) 2 


6 


(3 6 ) 


(l + v) 2 {l-4v + v 2 ) 2 



Table 2: The polynomials P\{v) as denned in equation ([7]), for each Archimedean lattice. 

for spin / A = F^/N for the lattice A, can finally be written as: 

f A /k B T = -p\og2- Z -\og{l+v 2 )- 1 - j -^^logdet [1 - vW A (k x , k y )} . (6) 

The matrices Wa(^, k y ) are m x m matrices with m = zp. Here z is the coordination 
number of the lattice A, while p is the number of lattice sites in the unit cell of A. The 
values of m, z and p for the Archimedean lattices A are reported in Table 1. 



4 Critical temperatures 

At the phase transition the free energy becomes singular, from equation ([6]) we see that 
this can happen only if the argument of the logarithm becomes zero. For each A this can 
happen only for k x = k y = 0. To every lattice we associate the polynomials: 

P A («) = det(l--uW A (0,0)) . (7) 

The critical values of the high temperature parameter < v c < 1 are the real solution of 
the equation 

Pa(v) = 0- (8) 

The polynomials P\{v) for the various Archimedean lattices are listed in Table 2. For all 
the lattices we considered, we found only one real solution of equation (JHJ) in the interval 
(0,1). In all, but in the cases of the lattices (3 4 ,6) and (3 2 ,4, 3,4), equation ([8]) has an 



6 



z 


A 


V c 


k B T c /J 


k B T*/J 


3 


(3, 12 2 ) 


1 _|_ 1 In _i_ 5^ 


1.2315 


5.0071 




(4,6,12) 


/ 5+3^- a/ 44+26^ 
V 2 


1.3898 


4.1363 




(4,8 2 ) 


1 + V 2 


1.4387 


3.9310 




(6 3 ) 


1 


1.5186 


3.6410 


4 


(3,4,6,4) 


1 , 

2 V 2 2 


2.1433 


2.4055 




(3,6,3,6) 


j / /TT /-- 

1 , / V3 i V3 

2 V 2 2 


2.1433 


2.4055 




(4 4 ) 


v 7 ^- 1 


2.2692 


2.2692 


5 


(3 4 , 6) 


0.344296 


2.7858 


1.8757 




(3 3 ,4 2 ) 


i 

3 


2.8854 


1.8205 




(3 2 ,4,3,4) 


0.32902 


2.9263 


1.7992 


6 


(3 6 ) 


2- >/3 


3.6410 


1.5186 



Table 3: For each Archimedean lattice, the exact critical value of the high temperature 
parameter v c , of the Curie temperature T c and of the dual Curie temperature T*, are 
given. The Curie temperatures T c and T* are calculated from v c using equations 
and ( froj) . The values for the lattices (4,6, 12), (3,4,6,4) and (3 4 , 6) were not yet known 
exactly and agree, within errors, with the numerical values of I2T] . The critical values 
for the lattices (3 3 ,4 2 ) and (3 2 ,4,3,4) agree with the exact results of [TS]. Note that the 
lattices (3,4,6,4) and (3,6,3,6) have the same critical temperature. 

explicit algebraic solution. The critical values v c are listed in the collum 3 of Table 3. 

The critical temperatures T c are related to the high temperature parameter v c by the 
following relation: 

Tc = J/kB (9) 
c arctanhf c 

The values we found are given in the collum 4 of Table 3. The exact Curie temperatures 
for all the lattices, except for (4,6,12), (3,4,6,4) and (3 4 ,6), were already known in the 
literature. In particular, they can be found in the following references: (4 4 ) in [21 [31 H], 
(3 6 ) and (6 3 ) in [6], (3, 12 2 ) in Q3], (4, 8 2 ) in Q3H Q3], (3, 6, 3, 6) in H2], (3 3 , 4 2 ) in 
[19] and (3 2 ,4,3,4) in P3 Q21 HI]. 

We have compared our exact results for the lattices (4,6,12), (3,4,6,4) and (3 4 ,6) 
with the numerical values available from Monte Carlo simulations [22| [21] and they agree 
within numerical errors. This is an independent confirmation of the correctness of our 
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analytical calculations for this new critical values. 

The critical temperatures for the lattices (3 2 ,4, 3,4) and (3 3 ,4 2 ) were calculated in 
[19J using the same combinatorial approach we are using, and agree with our findings. 
The fact that we used lattices where all edges are of unit length, contrary to [19J who 
used a different "embedding" of the lattices by constructing them from a square lattice 
adding diagonal links, shows that the critical temperatures depend only on the topology 
of the lattice and not on the particular representation of it. In fact, different embeddings 
of the same lattice lead to different random-walk matrices Wa(0, 0). This give rise to 
different polynomials Pa{v), that instead have the same roots in the interval (0, 1). The 
fact that the Curie temperatures depend only on the topology of the lattice is understood 
from the fact that the function $a(u), which counts the number of closed even graphs, is 
a topological property of the lattice which is independent of the embedding. 

Note that the lattices (3,4,6,4) and (3,6,3,6) have the same Curie temperature. 
This comes from the fact that the relative polynomials -P(3,4,6,4)(^) and P(3,e,3,6){ v ) have 
a common factor with root in the interval (0,1). This may suggest that the critical 
probabilities p c for correlated percolation (also know as bootstrap percolation) on this 
lattices are equal. 

Using the critical values for the high temperature parameter v c found for the Archimedean 
lattices, we can calculate the Curie temperatures for the Laves lattices from the relation: 

rp* _ _2J/kB_ 
\ogv c 

which follows from duality, as shown in |2] . The values of the Curie temperatures for the 
dual lattices A* are reported in collum 5 of Table 3. 



5 Conclusions 

We calculated the exact Curie temperatures for the two dimensional Ising model on all 
Archimedean and Laves lattices using the combinatorial approach of Feynman and Vod- 
vicenko. This results are summarized in Table 3. In the particular case of the lattices 
(4,6,12), (3,4,6,4), (3 4 ,6) and their duals [4,6,12], [3,4,6,4], [3 4 ,6] our results were 
known in the literature only approximately from Monte Carlo simulations. This numer- 
ical results agree, within errors, with our analytical findings. The universal character of 
the phase transition is reflected in the singularity of the free energy, which is logarithmic 
(a = 0) for all lattices. The Curie temperatures are not universal and we checked that 
they depend only on the topology of the lattice, but not on the particular embedding 
of it. The table of critical temperatures we compiled can be a useful reference to future 
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numerical and theoretical work on the Ising model on a large class of two dimensional 
lattices. 



Acknowledgments 

I would like to thank Prof. K. Binder for reading the manuscript and for useful discussion 
and suggestions. 

References 



[i] 

[2] 
[3 
[4 
[5 
[6 
[7 
[8 

[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 



R. Peierls, Proc. Camb. Philos. Soc. 32, (1936) 477. 

H. A. Kramers and G.H. Wannier, Phys. Rev. 60, (1941) 252. 
L. Onsager, Phys. Rev. 65 (1944) 117. 

B. Kaufman, Phys. Rev. 76 (1949) 1232. 

C. N. Yang, Phys. Rev. 85 (1952) 808. 
G.H. Wannier, Phys. Rev. 79 (1950) 357. 

K. Kano and S. Naya, Prog. Theor. Phys 10 2 (1953) 158-172. 

I. Syozi in Phase Trasitions and Critical Phenomena, C. Domb and M. S. Green eds. 
(1972) Wiley, New York, v. 1 pag. 269. 

T. Utiyama Prog. Theor. Phys. 6 (1951) 907. 

V. G. Vaks, A. I. Larkin and Yu. N. Ovchinnikov, JETP 22 820 (1966). 
Viktor Urumov, J. Phys. A: Math. Gen. 35 (2002) 7317-7321. 
J. Strecka, Phy. Lett. A 349 (2006) 505-508. 



V. Matveev and R. Shrock, J. Phys. A 28 (1995) 5235, arXiv: |hep-lat/ 9412105 
M. Kac and J. C. Ward, Phys. Rev. 88 (1952) 1332. 

R. P. Feynman, Statistical Mechanics, Benjamin/Cummings, Reading, MA, 1972 
N. V. Vdovichenko, Sov. Phys. JETP 20 (1965) 477-9. 
S. Sherman, J. Math. Phys. 1 202 (1960). 



9 



[18] T. Morita, J. Phys. A: Math. Gen. 19 (1986) 1197-1205. 
[19] C.J. Thompson and M.J. Wardrop, J. Phys. A 7, L65 (1974). 

[20] B. Griinbaum and G. Shepard, Tilings and Patterns: an Introduction, Freeman, 1989. 

[21] K. Malarz, M. Zaborek and B. Wrobel, TASK Quarterly 9 (2005) 475, arXiv: 
cond-mat/0506436[ 

[22] F. W. S. Lima, J Mostowicz and K. Malarz, arXiv: 1001.3078 [cond-mat.stat-mech]. 

[23] P. Fonseca and A. Zamolodchikov, J. Stat. Phys. 110 527 (2003). 

[24] G. Mussardo, Statistical Field Theory, Oxford University Press, 2009. 

[25] M. Kardar, Statistical Theory of Fields, Cambdrige University Press, 2007. 



10 



